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Abstract 

We study the propagation properties of the solutions of the finite-difference space semi-discrete wave equation on 
an uniform grid of the whole Euclidean space. We provide a construction of high frequency wave packets that 
propagate along the corresponding bi-characteristic rays of Geometric Optics with a group velocity arbitrarily 
close to zero. Our analysis is motivated by control theoretical issues. In particular, the continuous wave equation 
has the so-called observability property: for a sufficiently large time, the total energy of its solutions can be 
estimated in terms of the energy concentrated in the exterior of a compact set. This fails to be true, uniformly on 
the mesh-size parameter, for the semi-discrete schemes and the observability constant blows-up at an arbitrarily 
large polynomial order. Our contribution consists in providing a rigorous derivation of those wave packets and 
in analyzing their behavior near that ray, by taking into account the subtle added dispersive effects that the 
numerical scheme introduces. 

Resume 

On etudie les proprietes de propagation des solutions de l'equation des ondes semi-discretisee en espace par 
differences finies sur une grille uniforme dans tout l'espace euclidien. On realise une construction de paquets d'ondes 
concentres a haute frequence qui se propagent le long des rayons bicaracteristiques de l'Optique Geometrique a une 
vitesse de groupe arbitrairement petite. Notre analyse est motivee par la theorie du controle. Plus precisement, 
l'equation des ondes continue verifie la propriete d'observabilite : pour un temps suffisament grand, l'energie 
totale de ses solutions peut etre estimee en fonction de l'energie localisee a l'exterieur d'un ensemble compact. 
Cette propriete n'est pas verifiee de maniere uniforme par rapport au pas de discretisation pour le schema semi- 
discret pour un temps fini quelconque, si bien que la constante d'observabilite semi-discrete diverge avec un taux 
polynomial arbitraire. Nous donnons une construction precise de ces paquets d'ondes et decrivons l'effet dispersif 
rajoute que le schema numerique genere autour du rayon de propagation. 
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II est bien connu (cf. [11]) que les solutions du probleme de Cauchy pour l'equation des ondes continue 
(i-dimensionnelle verifient la propriete d'observabilite : pour un temps sumsament grand, l'energie totale 
des solutions peut etre estimee en fonction de l'energie localisee a l'exterieur d'un ensemble compact. Par 
la methode d'unicite de Hilbert (HUM) introduit dans [5], cette propriete est equivalente a un resultat 
de controlabilite exacte par un controle localise dans le complementaire du meme ensemble compact. 

Dans cet article, on considere le probleme de Cauchy associe a l'equation des ondes d-dimcnsionnelle 
semi-discretisee en espace par un schema centre en differences finies dans un maillage uniforme. Nous 
avons deux objectifs. Le premier est de montrer l'existence de paquets d'ondcs concentres autour d'un 
nombrc d'onde fixe a priori, a haute frequence, propageant avec une vitesse de groupe arbitrairement 
petite autour des rayons bicharacteristiques du schema semi-discret et done peu visiblcs depuis le domaine 
d'observation. On deduit dans [3] que la constante d'observabilite semi-discrete diverge au moins d'une 
maniere polynomialc arbitraire, et ce pour tout temps d'observabilite fini. Lc dcuxicmc objectif est de 
donner une forme asymptotique precise pour ces paquets d'ondes, de fagon a mettre en evidence la 
dispersion numcrique qui n'apparait pas dans le modele continu. 

Notre resultat complete la littcraturc cxistante sur l'observabilite et le controle des schemas d'approxi- 
mation de l'equation des ondes. En particulier, dans [10], pour la semi-discretisation de l'equation des 
ondes par differences finies dans des domaines bornes comme les cubes d-dimensionnels, il a ete prouve 
que l'incgalitc d'observabilite explose au sens ou la constante d'observabilite explosc lorsque le pas du 
maillage tend vers zero. Dans [7], on montre aussi que pour l'equation des ondes unidimcnsionnelle sur 
un intervalle borne observee par le bord, la constante d'observabilite diverge de maniere exponenticllc. 
Ce type de phenomene pathologique apparait aussi au niveau des estimations dispersives de Strichartz 
associees a l'equation de Schrodinger semi-discretisee en espace par differences finies, cf. [4]. 

L'avantage de notre construction de paquets d'ondes numeriques a haute frequence est qu'elle peut etre 
generalisee a d'autres approximations plus sophistiquees de l'equation des ondes comme celles obtenues 
par les methodes de Galerkin discontinus ou par des elements finis classiques d'ordre superieur, ou le 
symbole "numerique" en Fourier du Laplacien est une matrice donnant lieu a plusieurs relations de 
dispersion, cf. [6]. 

En outre, les paquets d'ondes que nous construisons sont aussi localises dans l'espace physique et 
peuvent done etre adaptes pour des problemes aux limites dans des domaines bornes avec des conditions 
aux limites diverses. 



1. Introduction and problem formulation 

For a finite time T > 0, consider the Cauchy problem associated to the d-dimensional wave equation: 



The problem (1) is well posed in fl' 1 (R d ) x L 2 (R d ) (H 1 ^ 4 ) is the completion of C^(R d ) with respect 
to the semi-norm || • d = j|V • H^RdO- The total energy of its solutions is conserved in time and is 




(1) 



defined by E(^,^) = (\m)\\l. IRd , + || W)|| 2 L2f W/2. 
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For all finite time T > 2, there exists a constant C(T) > such that, for all (<£°, (j) 1 ) £ H x {M. d ) x L 2 
the following observability inequality holds (see [11]): 

T 

E(<j>°, (p 1 ) ^ C(T) f ( \d t 4>{x,t)\ 2 dxdt. (2) 



o nd 

The observation region D, d := M. d \ B d (0,l) is the complement of the d-dimensional unit ball. This 
observability problem is motivated by controllability issues, since, by means of the Hilbcrt Uniqueness 
Method (HUM) introduced in [5], the inequality (2) is equivalent to the fact that, for all T > 2 and all 
initial data (u°, u 1 ) e fl' 1 (R d ) x L 2 (R d ), there exists a control function / <= L 2 (Vt d x (0, T)) such that the 
solution of the inhomogeneous Cauchy problem 

d 2 u -Au = f Xa *{x), x eK d ,te (0,T] 

u(a:,0) = it (a:), d t u(x,0) = u x {x), x £ R d 

satisfies u(x,T) — d t u(x,T) = for all x £ R d , where Xn d is the characteristic function of the set tt d . 

These issues are by now well understood for the continuous wave equation. In particular, observability 
inequalities of the form (2) and several variants hold under the sharp Geometric Control Condition 
(GCC) (cf. [1]) requiring that all rays of Geometric Optics enter the observation subdomain during 
the observability time. As shown in [8], when the GCC is not satisfied, the observability property fails 
because of the existence of Gaussian beam solutions localized around a bi-characteristic ray that escape 
the observation region during the time interval [0, T]. 

This paper is devoted to analyze these issues for the finite difference scmi-discrctc schemes. Given a 
mesh size h > 0, we define an uniform grid of the whole Euclidean space by Xj = h], j £ 7L d . We also 
introduce two discrete operators: the gradient = (d£ k )k=i,— ,d and the Laplacian Ah, with 



1=1 

where / = (/j)j e z<i is any sequence and (e;)f =1 is the canonical basis in R d . 

We consider the finite difference space semi-discretization of the wave equation (1): 

d 2 rf(t) - A fc #(t) =0, j £ Z d ,t e (0,T] 

4(0) = 4t >0 , $#(<)) = j e z". 

The problem (5) is well-posed in k 1 x £ 2 , with 



(4) 



(5) 



£ 2 = {f s.t. \\f\\l := h d \hf < ^} ^d ft 1 = {f s.t. \\~fwl := ||V+?|| £ 2 2 = £ \\d+J\\% < 



fe=i 



The total energy associated to its solution is conserved in time, being defined by 

E h ($ h ' ,^) = (0 l (t)\\ 2 hl + \\dJ h (t)\\ 2 2 )/2. 
For a fixed T > 0, consider the semi-discrete version of the observability inequality (2): 

W h,0 >"? M )<c , /iCn [ T \\d t ^(t)\\% m dt, with \\d t ^(t)\\% m = h d m(t)\ 2 - w 

For all finite T > and mesh size h > 0, it is easy to see that (6) holds for Ch(T) sufficiently large. 
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For the finite difference and Pi-classical finite element semi-discretizations of the wave equation on 
particular bounded domains like d-dimensional cubes, it is well-known from [10] that the corresponding 
observability constant blows-up as h — > 0, because of the pathological behavior of the spurious high 
frequency numerical solutions. In the context of the finite difference semi-discretization of the 1 — d 
wave equation on a bounded interval, the observability constant is known to blow-up exponentially as 
/i40 for all T > (cf. [7]). As shown in [4], similar pathological high frequency phenomena have been 
observed in the context of the Strichartz dispersive estimates for the finite difference approximations of 
the Schrodinger equation. 

In view of the formulation of the observability problem under consideration, our goal here is twofold. 
Firstly, to show the existence of numerical waves that are concentrated in space-time and then to give its 
precise asymptotic form in order to illustrate the added effect of the numerical dispersion with respect to 
the classical continuous wave equation. These solutions propagate according to the group velocity (cf. [9] ) , 
a notion that our construction contributes to make it precise. The group velocity of the high frequency 
numerical waves may be arbitrarily small, as the analysis of the dispersion diagram below shows. 

Set 11^ := [— n/h, n/h] d . By applying the semi-discrete Fourier transform (SDFT) (see [9]) on (5) and 

denoting by <p h (^,t) the SDFT of the solution (j) h {t) of (5), one obtains the following second-order ODE 
depending on the parameter £: 

^,o) = ^°(0, a t ^,o) = ^ 1 (0,^e< 

where u)d,h{C) is the multi-dimensional dispersion relation associated to (5), ^h(0 = W Sfe=i sm2 (%^)> 

£ = (&)i<fc<d e ii d h - 

The semi-discrete rays of Geometric Optics corresponding to the semi-discrete problem (5) are straight 
lines of the form x^(t) = x ± Vw^/,^)*, x G R d , £ G 11^, propagating with a group velocity, defined as 
the gradient Vu>d,h(0- Note that | Vujd,h(0\ vanishes for all £ G {±w/h, 0} d \ {0}. This is in contrast with 
the behavior in the continuous case (1), where the dispersion relation is cjd(£) ~ |£| and the velocity of 
propagation along all rays is IVw,^)! = 1, for all £ G R d . 



2. Statement of the main result 

Let T > be given. Choose an arbitrary x* G B d (0, 1) and consider a wave number 770 = fr£o S n^\{0} 
such that the corresponding semi-discrete ray starting at x* and traveling with velocity |Vw<2,/i(£o)l — 
IVw^i^o)! does not enter the observation region in time T, i.e. \x* — iVixVa^o)! < 1, f° r all t G [0,T]. 

Consider 7 := 7(^1) > such that 

7 >> 1 and hj « 1. (8) 
For G <S(M d ), the class of Schwartz functions in R d , </> being its continuous Fourier transform, consider 
the semi-discrete wave equation (5) with initial data ( <j) h '° 7 4> h ' 1 ) whose SDFT is given by 

M (O - -^TT-^r^ £ f — ) cxpHx* • K - 6>))Xnj[(0 and ^(0 - iu, d>h (Q^°(Q. (9) 
iUd,h(£) V 7 V 7 / h 

These data, in the physical space, correspond, roughly, to highly oscillating Gaussian profiles. As shown 
by S. Ervedoza in [3], by a stationary phase like argument, one can show that for all /3 G K + , there exists 
a constant Cp((j>,T) such that the observability constant C^(T) in (6) satisfies Ch{T) ^ Cp{4>,T)"fP . 

The initial data (9) are concentrated high frequency wave packets at the wave number £ = £ = Vo/h 
with a width of order I/7, with 7<<l//i^ooas/i^0. This limits drastically the spread of the wave 
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packet and allows to concentrate its energy around the chosen ray and to reduce the added dispersive 
effects that the non-trivial Hessian matrix of the dispersion relation introduces. Simultaneously, by the 
uncertainty principle, this forces the wave packet, in the physical space, to have a spread factor of order 
I/7 >> h, which is asymptotically larger than the characteristic size h of the mesh. This is natural, in the 
sense that the numerical effects are only detected when an infinite number of nodes enter asymptotically 
in the determination of the data of the solution. 

This example of wave packet clearly illustrates the classical effect due to the group velocity that is 
plotted in Figure 1 below. We plot in red the initial velocity ft(x) = exp(— "f\x\ 2 /2) exp(z£ £) (we also 
take i\$\ft = ft), with h = 0.005, £ /i = 19tt/20, 7 = h~ a , 2a = 0.75 in d = 1 (left) and d = 2 (right), 
which coincides, up to an exponentially small error (with respect to h), with the discrete initial data (9). 
In blue, we plot the time derivative of the solution of the continuous wave equation (1) and in black the 
one corresponding to the semi-discrete wave equation (5) with initial data (9) at time t — 1. For d = 2, we 
represent only the projection of these (continuous and discrete) solutions on the x-plane. This experiment 
shows that, as theory predicts, the semi-discrete wave packets propagate with a velocity that is much 
smaller than the one of the continuous wave equation. 






Figure 1. Transport of continuous versus discrete high frequency wave packets in dimension d = 1 and d = 2. 

In the following, we analyze more precisely the behavior of these wave packets. To do this, we appro- 
ximate the analytic dispersion relation uJd,h{£.) by a polynomial one. 

Set IE — £0 : = {£, — £.0; £ € n^}. The time derivative of the solution of (5) with initial data (9) is given 
by the following wave packet 

-d 



1 



(2tt) c 



0(f/7)Xn£-£ o (O ex P( itw d,'i(£ + £o))exp(^- (x } -x*))d(. 



(10) 



To simplify the presentation and, without loss of generality, we set £0 = and x* — 0. Since (j> € S(M. ), 
we may also neglect the characteristic function in (10). Setting w(£) = w<j i(£ + 770) and x — — x*, the 
wave packet (10) can be written as 



u(x, t) 



(27T) 



1 7 



) exp ^— u)(£h)j exp(z£ • x) d£. 



(11) 



In view of the analyticity of u>(r}), we can split it as oj(t)) = L(rj) + D(rj) + R(ij), where L(ry) = 
w(0) + Vw(0) • rj is the linear part of w, D is the second order term in the Taylor expansion of w about 
rj = and R is the corresponding reminder, given explicitly by 

1 

E ~D a u(0)r 1 a ,R(r,)= ]T ~ V a J (1 - X) 2 D a u(\ V ) d\. 



H=2 



H=3 



Factoring out the time dependent complex exponential generated by the zero order term in the Taylor 
expansion of the dispersion relation, we may decompose u as uexp(— itLO^\{Q)/h) = v + v R , where 

U(M)= (2^F / ^) cx p( l l D ^ h )) c M^-(x + tVoj(0)))d^ 



and 



R d 

Our main result determines the conditions on 7 and <j>, such that is a reminder term. 
Theorem 2.1 For all t > and </> € 5(IR <i ) 7 the following estimate holds: 



+M12 



L2(Rd) < C$)h 4 7 H 2 , (12) 



wfeere C(0) = C||| • l^ll^^/IMI^) and C = ( ^ sill^^H^WOM))) 

For a finite time interval £ € [0, T], one can guarantee that v R is small in the sense of (12) when 

7 /i 2/3 << 1, (13) 

which is a more restrictive condition on 7 than (8), guaranteing that the energy concentrated outside 
the ray is polynomially small. Indeed, (12) yields an asymptotic description of the solution globally in 
space-time, i.e. u ~ exp(itco(0)/h)v. Intuitively, the scale (13) is motivated by the fact that \R(£h)/h\ < 
VCh 2 \£\ 3 . For this to be asymptotically small, the width of the Fourier transform of the profile has to be 
limited by (13). 

The function v is a solution of the following PDE: 

d t v = Vw(0) • V x v -ih^2 D a oj(0)D%v. (14) 

M=2 

This is a transport equation perturbed by an asymptotically small (as h tends to zero) Schrodingcr 
like second order term. The following result emphasizes that there exists solutions of (14) for which the 
relevant scale is 

7/i 1 / 2 = 1. 

Theorem 2.2 There exist solutions of (14) admitting the following asymptotic expansion: 



v(x, t) 



t^fai&M*- < 15 > 

j=0 v 7 

where 770 G K d is a fixed wave number, cj 2 (?7) = Vw(0) ■ + , =2 ^D a uj{Q)rj a is the dispersion relation 
corresponding to (14) and (aj(x,t))j e ^ solve the following system of PDEs: 

d t a = Vw(0) • V x a Q and 8 t a j+1 = Vw(0) • V x a J+1 - 1 £ ^D a w(0) i^ofD^aj^j € N. 

\a\=2 ' 0</3<a 

By taking the support of the initial datum of ao to be compact, one can observe that v is concentrated 
along a neighborhood of the ray of width \fh. This scale \[h is critical due to the added dispersion 
that the Schrodinger like term introduces. Once the asymptotic expansion of v is given as in (15), one 



immediately gets that of u and therefore of d t 4> and of </> . In this way, we get the asymptotic form of 
the high frequency wave packet. 

This kind of expansion can be further developed, incorporating higher order terms of the Taylor ex- 
pansion of the dispersion relation and a multiple-scale ansatz. This issue will be developed in detail in 
[2]- 
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